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The reaction v e + p — > e + + n is very important for low-energy (Ev < 60 MeV) antineutrino 
experiments. In this paper we calculate the positron angular distribution, which at low energies is 
slightly backward. We show that weak magnetism and recoil corrections have a large effect on the 
angular distribution, making it isotropic at about 15 MeV and slightly forward at higher energies. 
We also show that the behavior of the cross section and the angular distribution can be well- 
understood analytically for E v < 60 MeV by calculating to 0(1/M), where M is the nucleon mass. 
The correct angular distribution is useful for separating u e +p — > e + + n events from other reactions 
and detector backgrounds, as well as for possible localization of the source (e.g., a supernova) 
direction. We comment on how similar corrections appear for the lepton angular distributions in 
the deuteron breakup reactions v e + d — > e + + n + n and v e + d — > e~ +p + p. Finally, in the reaction 
Ve + P — * e + + n, the angular distribution of the outgoing neutrons is strongly forward-peaked, 
leading to a measurable separation in positron and neutron detection points, also potentially useful 
for rejecting backgrounds or locating the source direction. 
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I. INTRODUCTION 

Inverse neutron beta decay, v e + p — > e + + n, is the 
reaction of choice for the detection of reactor antineutri- 
nos, crucial to neutrino oscillation searches. It is also, by 
far, the reaction giving the largest yield for the detection 
of supernova neutrinos. The LSND and KARMEN @] 
experiments use antineutrinos from /x + decay at rest 
to search for the oscillation appearance of v e events, de- 
tected by this reaction. There are also searches for v e 
antineutrinos from the Sun. 

In many of these applications, in particular those based 
on detection by Cerenkov radiation, one can determine 
the direction of the outgoing positron. It is therefore of 
interest to consider the angular correlation between the 
incoming antineutrino and outgoing positron directions, 
and its energy dependence for E v < 60 MeV, relevant for 
the above studies. If the source direction is known, the 
angular correlation can be used to help separate these 
events from other reactions or detector backgrounds. If 
the source direction is unknown, as possibly for a Galactic 
supernova, then the observed angular distribution may 
help to locate the source. 

For low antineutrino energies, the positron angular dis- 
tribution is well-described by the form 



dg 
dcosf? 



1 + v e a(E v ) cosO , 



(1) 



where 6 is the angle between the antineutrino and 
positron directions in the laboratory (where the proton 
target is assumed to be at rest) and v e is the positron 
velocity in c = 1 units. At higher energies, terms pro- 
portional to higher powers of cos 9 appear, and at the 
highest energies, all reaction products are strongly for- 
ward simply by kinematics. It is convenient to describe 



the angular distribution by the average cosine, weighted 
by the differential cross section,since that is always well- 
defined. In the limit that Eq. (|l|) holds, 



-v e a(E u ) , 



(2) 



where the E^-dependence of (cos 9) is suppressed in the 
notation here and below. Except near threshold, where 
v e becomes very small (but nonzero in the lab), the v e 
factor is nearly unity and can be ignored. 

In the limit where the nucleon mass is taken to be 
infinite, i.e., zeroth order in 1/M, the asymmetry coeffi- 
cient a is independent of E v and would be the same for 
Ve + V ~~ * e+ + n an d v e + n — ► er + p. (Since there are 
no free neutron targets, the latter cannot be directly ob- 
served. We return to this point below in discussing neu- 
trino and antineutrino reactions with deuterons.) Then 
is given simply by the competition of the non-spin- 
flip (Fermi) and spin-flip (Gamow- Teller) contributions, 
and is 



,(°) 
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(3) 



and thus the angular distribution of the positrons is 
weakly backward. We have defined the vector and axial- 
vector coupling constants by / = 1, g = 1.26. 

In the following we will consider how (cos 9) is modified 
when weak magnetism and recoil corrections of (9 (1/M) 
are kept. The effect of these terms on the total cross sec- 
tion was calculated in Refs. [|]]|] (see also Ref. Ej), where 
it was found that they, in particular the weak magnetism, 
reduce the cross section by a noticeable amount. In this 
paper, we will show that the positron angular distribu- 
tion is changed even more, including the sign of (cos 9). 
The effect is so large in part because of the accidental 
near-cancelation in . 
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The general form of the differential cross section, valid 
to all orders in 1/M but neglecting the threshold effects 
(and hence only valid for energies far above threshold) , is 
well-known ||. For the relevant energies E v < 60 MeV, 
it is instructive and sufficient to consider in detail just the 
terms of first order in 1/M. Here and below, 1/M will be 
taken to refer to all terms of that form, with E v /M being 
dominant among them. Moreover, using these results, we 
show how to extend the formula of Ref. || to low energies, 
so that it merges smoothly with the correct expression 
near threshold. 



II. THE POSITRON ANGULAR DISTRIBUTION 

A. Differential cross section: expansion in powers of 

1/M 

We begin with the matrix element of the form 



M 



G f cos 8c 

X [VuJ^{l ~ 75) V, 



7/J - 7^753 - 



2M 



(4) 



where / and g are given above, the anomalous nu- 
cleon isovector magnetic moment is defined with f% — 
Mp — Mn = 3.706, and cos 9c = 0.974. In the most gen- 
eral case, the coupling constants are replaced with form 
factors that vary with q 2 ; we neglect this variation as it 
is 0(E 2 /M 2 ). The four-momentum transfer q 2 is related 
to the laboratory scattering angle 9, which in turn is re- 
lated to the outgoing positron energy E e (again in the 
laboratory) by the relations 



q 2 = m 2 - 2E v E e (l - v e cos8) 
= {M 2 n - M 2 ) - 2M p (E v - E e ) 



(5) 
(0) 



Some other useful kinematic relations are given in the 
Appendix. We can now use the standard rules and evalu- 
ate the differential cross section accurate to a given order 
in 1/M. 

At zeroth order in the positron energy is 



E^ = E u - A, 



(7) 



where A = M n — M p . At each order in 1/M, we define the 
positron momentum p e = \J E 2 — m 2 and the velocity 
v e = p e /E e . The differential cross section at this order is 



Vdcos(9 J 

CO 
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(/ 2 + 3. 9 2 ) + (/ 2 - ff >fcos0 E?V?K (8) 



The normalizing constant <tq, including the energy- 
independent inner radiative corrections, is 



cr = 



G 2 F cos 2 9 C ^ 



inner / J 



(9) 



where A- 

v inner 



0.024 |7j. This gives the standard result 
for the total cross section, 
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The energy-independent inner radiative corrections af- 
fect the neutron beta decay rate in the same way, and 
hence the total cross section can also be written 



(0) _ 2n 2 /ml (0) (0) 

"tot — j-p "e fe 
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where r n is the measured neutron lifetime and fp S — 
1.7152 is the phase space factor, including the Coulomb, 
weak magnetism, recoil, and outer radiative corrections, 
but not the inner radiative corrections M . The cross sec- 
tion normalization was measured in Ref. H and found to 
be in agreement with the expectation from the neutron 
lifetime at the 3% level. The (small) energy-dependent 
outer radiative corrections to a tot are given in Refs. ||Q. 
The outer radiative corrections to (cos 9) should largely 
cancel in the ratio of the cross section weighted with cos 9 
to the cross section itself, and so are not considered fur- 
ther here. 

At first order in 1/M, the positron energy depends 
upon the scattering angle and is 



E^ = 2#> 



l-^ ( l- v (o) cos( 



where y 2 = (A 2 — m 2 )/2. In factors of the form 1/M, we 
use the average nucleon mass; using 1/M versus 1/M p 
leads to an ignorable difference of C(l/M 2 ). The differ- 
ential cross section at this order (after a lot of tedious 
algebra) is 
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where 

r = 2(f + f 2 )g 

+ (f+9 2 ) 
+ (f 2 +3g 2 ) 

+ (f-9 2 ) 
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For the dominant term (the first square brackets in 
Eq. (HH)), the cosine-dependence of E e , p e , and v e must 
be taken into account at 0(1/M), while for the subdom- 
inant terms, they may be taken as functions of E v alone. 
Note that the terms above, which can be large 

near threshold, are canceled by p^ in the phase space 
factor. Our result for the total cross section, calculated 
from Eq. (pd]), supercedes the result of Ref. which 
did not include the 0(1/M) corrections to the Jacobian 
dt/d cos 9. 

Unless the electron mass is negligible (see below), it is 
not convenient to analytically expand Ee , p^\ and 
in powers of 1/M, and instead we evaluate the total cross 
section a^ t numerically. In the upper panels of Figs. 1 
and 2 we show the total cross section versus E v . We di- 
vided the plots into two energy regimes: Fig. 1 to show 
the threshold region, relevant for, e.g., reactor experi- 
ments; and Fig. 2 to show the global behavior relevant 
for, e.g., supernova or muon decay at rest experiments. 
The solid line is the result at first order in 1/M, given by 
numerical integration of Eq. (114) . The short-dashed line 
is the result at zeroth order in 1/M, given by Eq. ( |TT| ) . As 
expected, these results agree at the very lowest energies. 
However, with increasing energy the 1/M corrections be- 
come more and more important, reducing the total cross 
section. 
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FIG. 1. Upper panel: total cross section for u e +p — > e + +n; 
bottom panel: {cos 8} for the same reaction; both as a func- 
tion of the antineutrino energy. The solid line is our 0(1/ M) 
result and the short-dashed line is the 0(1) result. The 
long-dashed line is the result of Eq.(3.18) of Ref. ||, and the 
dot-dashed line contains our threshold modifications to the 
same. The solid and dot-dashed lines are not distinguishable 
in this figure. The inner radiative corrections are included 
(see the text), but the outer radiative corrections are not (see 
Refs. §,§). 



B. Differential cross section: the high-energy limit 

Far above threshold, our result can be compared to 
Eq. (3.18) of Ref. ||; as noted, that formula neglects 
A but contains all orders in 1/M. At low energies, the 
neglect of the threshold is a large effect, as shown by the 
long-dashed line in the upper panels of Figs. 1 and 2. We 
have modified Eq. (3.18) of Ref. j(| to take into account 
the largest contributions of the threshold effects. First, 
the exact kinematics (see our Appendix), including A, 
should be used to evaluate q 2 and s — u in that formula. 
With no further modification, that formula does not have 
the correct low energy limit (determined by comparing to 
the results above). By direct comparison, the only other 
A-dependent correction to the formula of Ref. |(| at 0(1) 
is the modification 

C{q 2 ){s - u) 2 C(q 2 ){s - u) 2 - C(q 2 ) AM 2 A 2 . (16) 

As shown by the dot-dashed line in the upper panels 
of Figs. 1 and 2, this corrects the result of Ref. [|| so 
that it takes into account the threshold. There may be 
additional corrections of order A/M necessary, but by 
the numerical results, they are evidently small. 

We took the form factors in Eq. (||) to be constants, 
which is reasonable for the energies considered. At 
higher energies, the form-factor variation with momen- 
tum transfer must be properly included, as done in 
Note that the form- factor variation in Ref. H is incorrect, 



as it attributes a dipole behavior to the coefficients which 
appear directly in Eq.(^), and not to the momentum- 
transfer dependent linear combinations of them known 
as the Sachs form factors; see Ref. |J. 

The plotted results show that our a^H agrees well with 
the modified results of Ref. || , particularly at low ener- 
gies. While lowest order in 1/M is clearly not enough, 
the first order in 1/M is, justifying our neglect of higher 
orders. One can see that the short-dashed line differs 
from the others substantially already at E v > 30 MeV. 
This suggests that in the expansion in 1 /M the numerical 
coefficient multiplying the dominant term E u /M is quite 
large (ci —7). One can see that explicitly by examining 
Eq. @) below. 

C. Angular distribution 

At zeroth order, recall that 

(cosfl) (0) = t£V 0) /3 ~ -0.034u^ , (17) 

so the angular distribution is slightly backward, inde- 
pendent of energy (above the threshold region). From 
Eq. (|TJ), it is evident that the angular distribution will 
be modified by the 1/M corrections. However, the cor- 
rections in Ee , pi , and are not explicitly shown. 
At lowest order, the m e -dependent effects appear always 
as m 2 /E 2 , and for E v ~>b MeV, these may be neglected 
so that 
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where 

r = 2(/ + / 2 )J(2£:( 



A)(l - cos< 
+ (f+g 2 )[A(l + cos9)} 
+ (/ 2 + 3.g 2 ) [s(E^ + A)(l - cosfl) - A 

+ (/ 2 - g 2 ) [3(^ 0) + A)(l - cos 6<) - A] cos 6», (19) 

where in the latter two square brackets, we have also 
neglected terms +A 2 /Ei°\ Note that the (/ 2 + 3.g 2 ) and 
(f 2 — g 2 ) terms are modified from Eq. ( fTi]) by terms from 

the expansion of the phase space factor E^p^K 

It is now trivial to integrate over cos 9 and to determine 
a tot and the integral weighted with cos 9, which we call 
(da cos 9) tot ■ These can be written as 
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and 
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The coefficients a, (3, 7 can be immediately read off of 
Eq. (|l8|), since the COS6 1 integration is trivial. In order 
to continue working consistently to order l/M, we divide 
(da cos 9)tot/&tot analytically. Since ^(\/ot\ — — 7 is large, 
numerical division would improperly introduce higher- 
order terms. Then 
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It can be shown both analytically and numerically that 
by far the largest m e -dependent effect can be restored by 
the insertion of the term vi°^ as above, since the phase- 
space factors E e p e cancel in the definition of (cos 9) . This 
formula is an excellent approximation for (cos 9) from 
threshold even to E v ~ 150 MeV (though at that energy 
neither of Eqs. (|2(]) and (|2l] ) is individually valid, and the 
angular distribution is no longer of the form 1 + a cos 9). 
Keeping only the largest terms, we can also write 



4(/ + / 2 )g \ K 

(p + m) M 



-0.034t4°) + 2.4^ 
M 



(23) 
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FIG. 2. The same as Fig. 1, but over a larger range of 
antineutrino energy. The long-dashed and dot-dashed lines 
are nearly indistinguishable in the lower panel. 



The standard M — > 00 result is very small, (cos6 ) )(°) ~ 

— 0.034i4°' > - The large corrections depending on E v jM 
can be classified as being due to weak magnetism (de- 
pending on (/ + f-i)) and pure recoil (independent of the 
couplings). For v e + p — > e + + n, these add. For the 
reaction v e + n — > e~ +p, the sign of (/ + fa) is reversed, 
and the recoil and weak magnetism terms would nearly 
cancel. 

Our results for (cos 9) are shown in the lower panels of 
Figs. 1 and 2. Our main result at order l/M, given by 
Eq. (p2|), is shown as the solid line. The zeroth order re- 
sult is shown as the short-dashed line. The long-dashed 
line shows the result of Eq. (3.18) of Ref. [gj, which 
assumes A = 0. This is obviously poor in the thresh- 
old region. The dot-dashed line shows our modification 
of that formula to account for the largest A-dependcnt 
effects. Note that the upward "hook" at low energies 
(E v < 5 MeV) is caused by the finite electron mass. At 
those energies, v e < 1, and the average (cos 9) decreases, 
nearly vanishing as the antineutrino energy approaches 
its threshold value. 

These results for (cos 9) agree qualitatively with the 
earlier numerical results of Ref. |l(]], which noted that 
the (cos 9) vanishes near E v — 20 MeV, and becomes 
slightly positive at larger antineutrino energies. At lower 
energies, the results of Ref. [|l0| are inaccurate, presum- 
ably due to using the formula of Ref. H without the 
threshold modifications given above. 

Terms of the first order in l/M radically change 
(cos 9), including its sign. At high energies, the missing 
l/M 2 terms become important for a to t and (not shown) 
(da cos 9)tot- However, note that l/M 2 effects are negli- 
gible for (cos 9). Our result, given by the solid line in the 
lower panels of Figs. 1 and 2, contains no terms of order 
l/M 2 or higher. The modified result (with our correc- 
tions for A > 0) of Ref. Q, given by the dot-dashed line, 
contains all terms of order l/M 2 and higher. The agree- 
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ment is excellent, and both are approximately linear in 
E„/M. That is, for (cos0), there is a large cancelation 
of the higher order corrections. 



D. Charged-current deuteron breakup reactions 

Since there are no free neutron targets, the reaction 
v e + n — > e~ + p cannot be observed directly. However, 
since the deuteron is so weakly bound, we can at least 
qualitatively apply the weak magnetism and recoil effects 
calculated above to the reactions v e + d — > e + + n + n and 
v e + d — > e~ + p + p. For the considered energies, these 
reactions are pure Gamow- Teller transitions, and so the 
asymmetry is = —1/3. In both reactions, (cosd) 
will be made more positive by pure recoil corrections. To 
those, the weak magnetism correction adds for v e + d and 
subtracts for v e + d. 

In Fig. 3 we show the (cos 8) results calculated from 
the double-differential cross sections of Kubodera JjT| . 
His results are based on a complete calculation, includ- 
ing treatment of the deuteron wave function and meson- 
exchange effects 1 12 . The corrections due to the finite nu- 
cleon mass are evident. One can see that the two curves 
are not symmetric with respect to the M — > oo value 
(cos 9) — —1/9. As expected, the weak magnetism and 



recoil corrections act in the same sense for P e + d and the 
opposite sense for v e + d. 

The weak magnetism contribution can be analytically 
estimated from the amplitude squared in Ref. [fL3| . We 
estimate the recoil contribution so that the combined re- 
sult 



9 3 V 9g J M 



(24) 



provides a reasonable fit to the full numerical results in 
Fig. 3. 



III. APPLICATIONS OF THE POSITRON 
ANGULAR DISTRIBUTION 

A. SN1987A events 

Supernova 1987A was observed in two water- Cerenkov 
detectors, Kamiokande II and 1MB Q, with 12 and 
8 events, respectively. These events were presumably en- 
tirely due to v e + p —* e + + n, with an angular distri- 
bution of the form 1 + a cos 8. A well-known peculiarity 
of the SN1987A data is that the angular distributions of 
the detected positrons are apparently too forward, with 
(cos 8) = 0.34 in Kamiokande II and (cos 8) = 0.48 in 
1MB. Using the results of Fig. 2, evaluated at the ob- 
served average energies, and the correction for the 1MB 
angular bias |15|, we would expect only (cos 9) ~ in 
Kamiokande II and (cos0) ~ 0.08 in 1MB. 

The error on the mean (cos 8) is 
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FIG. 3. The average lepton cosine for the charged-current 
deuteron reactions, versus the neutrino or antineutrino en- 
ergy, using Kubodera's calculations. We plot only from 1 
MeV above threshold, so the w e -dependence at low energies is 
not shown. The small jitters are due to the coarse integration 
grid. 



6 ((cob 6)) = 



or 



(25) 



since \a\ <C 1. Thus the expected error on (cos 8) from 
just statistics is 0.17 for Kamiokande II and 0.20 for 1MB. 
The range of antineutrino energies contributing only neg- 
ligibly increases the error. Thus the experimental results 
for both Kamiokande II and 1MB deviate by +2-cr from 
the expectations. (The disagreement between the exper- 
imental results and the expectations is also discussed in 
Refs. JlClplI ), At the same time, however, after correct- 
ing for the energy difference and the 1MB angular bias, 
the two means are in good agreement with each other. 

It is generally assumed in the literature that the 
angular distributions of Kamiokande II and 1MB are 
consistent. For example, Ref. jl7j claims an 81% 
Kolmogorov-Smirnov probability that the distributions 
are the same. That test is primarily sensitive to differ- 
ences in means ]l8| ] , and so this confirms the agreement of 
the (cos 8) values above. However, the angular distribu- 
tion is also characterized by its variance (cos 2 8) — (cos 8) 2 , 



with the expectation being : 
ance is lfl9| : 



1/3. The error on the vari- 



5(<cos 2 0)-(cos0) 2 )~^= A /^-i 



0.30 



TV 



(26) 



The experimental result for the variance in Kamiokande 
II is 0.32, with expected error 0.09, and hence in excellent 
agreement with expectation. The experimental result for 
1MB is 0.11, with expected error 0.11, and hence a — 2-er 
deviation from expectation and, more importantly, from 
the Kamiokande II result. Thus, contrary to general be- 
lief, the Kamiokande II and 1MB angular distributions, 
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characterized here by their first two moments, are not 
consistent at the 2-cr level. 

It is possible that some of the observed forward events 
were due to neutrino-electron scattering, though the ex- 
pectation is only ~ 0.3 events for Kamiokande II and 
~ 0.1 events in 1MB (using the same supernova parame- 
ters as in Rcf. [^(J and the detector properties taken from 
Refs. |]l4| , |l5l| ) . Most other authors have also obtained 
an expectation of < 1 event per detector. Allowing n 
neutrino-electron scattering events out of a total of N 
events will change the expectations for the mean (cos 9) 
(increased by n/N) and the variance (cos 2 9) — (cos6*) 2 
(increased by 2/3(n/N) - {n/N) 2 ). Thus for N ~ 10 and 
the possible n = 1, the means would be somewhat im- 
proved (now each a +1.5-a deviation). The Kamiokande 
II variance would still be in agreement with expectation, 
though the 1MB variance would then be a — 3-tr devia- 
tion. 

As a general caution about the small- number statistics, 
one can consider, for the purpose of illustration, the effect 
of assuming that one backward event was missed. That 
is, to each data set we add one fake backward event. 
For Kamiokande II, the effect on both the mean and the 
variance is modest, but for 1MB there is a large effect, 
making the mean only a +l-a deviation and putting the 
variance in agreement with theory. Thus the statistical 
significance can be very sensitive to small fluctuations, so 
that the number of sigmas of deviation and the implicit 
confidence levels should be taken with caution. 

In conclusion, the Kamiokande II data seem to re- 
quire a 2-a statistical fluctuation in the mean, and the 
1MB data separate 2-cr fluctuations in the mean and 
the variance. It is difficult to explain the observed an- 
gular distributions, even taking into account the cor- 
rections of Fig. 2 and (somewhat implausibly) allowing 
~ 1 neutrino-electron scattering event in each detector. 
Given the inconsistency between the Kamiokande II and 
1MB angular distributions, it is probably not legitimate 
to combine them, thus weakening the argument for new 
supernova or particle physics that could have affected 
the angular distributions, as invoked in Refs. fl7| , |2l| , p2"f . 
While we have not explained the angular distributions, 
we have explicitly shown the perils of the small-number 
statistics and an apparent additional problem with the 
1MB results. 



B. Supernova antineutrinos 

A strong P e signal is expected in SuperKamiokande |25| l 
and other underground detectors from a future Galac- 
tic supernova (for the expected count rates see, e.g., 
Refs. [ p0|p4| ]). Is it possible to use the observed angu- 
lar distribution of the positron events to locate the di- 
rection of the supernova p5|? If the M — > oo limit were 
appropriate (i.e., Eq. (p~7|)) , then the positrons would be 
dominantly moving in the backward direction. However, 
the l/M corrections calculated above are very important. 



Folding in the expected Fermi-Dirac distribution of the 
incoming v ei and weighting cos 9 properly with the flux 
and cross section calculated here, we arrive at (cos#) ~ 
0.015, 0.025 and 0.034 for temperatures T = 4, 5, and 6 
MeV. Thus the positrons from supernova v e , with most 
probable temperature of about 5 MeV, will in fact be 
slightly forward and, moreover the asymmetry coefficient 
will sensitively depend on the antineutrino temperature, 
quite different from the naive expectation. For locating 
the supernova, the best strategy seems to be to concen- 
trate on the positrons of the highest energy. For T = 5 
MeV, about 25% of the signal will be above 30 MeV, and 
should have a noticeable forward asymmetry ((cos0) = 
0.056). Observation of the angular distribution of the 
higher energy positrons would constitute an important 
check of the supernova origin of the signal, and would al- 
low location of the supernova to about <5(cos 9) ~ 0.2 [EsJ . 

If the supernova direction is known, then knowledge 
of the positron angular distribution could be used to 
separate these events from other reactions. For exam- 
ple, if there are v r — > v e oscillations, then the reac- 
tion v e + 16 O — > e~ + 16 F can become important |26] ] 
(since the v T temperature is higher than the v e temper- 
ature and this reaction has a relatively high threshold). 
The outgoing electrons are somewhat backward. Note 
that the neutron in T> e + p — ► e + + n is not detected, 
and electrons and positrons are indistinguishable by their 
Cerenkov radiation. Therefore, the search for events from 
v e + 1& — > e~ + 16 F must be done statistically, by looking 
at the total angular distribution and looking for a back- 
ward excess over what is expected from v e + p — > e + + n 
alone. The calculation in Ref. [^7j used the naive positron 
asymmetry and would have to be revised. Since there 
would be fewer backward events than they expected, the 
sensitivity would be improved. 

In a heavy water detector like the Sudbury Neutrino 
Observatory |2£| , the angular distributions of the outgo- 
ing leptons in the reactions i> e -\- d — > e + + n + n and 
v e + d — > e~ + p + p could also be used to locate the su- 
pernova direction. Because of the low numbers of events, 
however, even with a naive asymmetry of = —1/3, 
the pointing resolution is only about S(cos9) ~ 0.5 p5[ . 
Taking into account the l/M effects weakens the positron 
asymmetry, and would degrade the pointing. However, 
the corrected angular distributions will still be quite im- 
portant for separating reactions. 

C. Search for solar antineutrinos 

Ref. [^9| discusses the possibility of searching in Su- 
perKamiokande for 9 e antineutrinos from the Sun, pre- 
sumably from v e — > v e oscillations, with E u ~ 10 MeV. 
The authors proposed that these events could be sepa- 
rated statistically by their angular distribution from the 
isotropic detector background and the forward-peaked 
solar neutrino events from neutrino-electron scattering. 
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However, from Fig. 2, the angular asymmetry is substan- 
tially weaker than they assumed, and in view of that, 
their derived limit would have to be modified. 



D. LSND results 

Another important application of the positron angular 
distribution is the search for neutrino oscillations by the 
LSND || and KARMEN § collaborations. The LSND 
collaboration reported evidence for — > v e oscillations 
following p + decay at rest. The evidence is based on the 
observation of 22 e + + neutron events with positron en- 
ergies between 36 and 60 MeV, presumably originating 
from the v e p interaction. The directions of individual 
candidate positron events have been measured with re- 
spect to the beam axis and the angular distribution is 
given in Fig. 21 of Ref. Q]. The measured (cos^) was 
found to be 0.20 ±0.13. 

From our Fig. 2, one would estimate that the expected 
value would be only about 0.08. However, there are 
important experimental corrections particular to LSND 
which must be taken into account. The physics basis of 
most of them is the simple fact that, for a fixed antineu- 
trino energy E v , the forward-going positrons have more 
energy than the backward-going ones. Thus, if there is 
a cut on positron energy, say E e + > 36 MeV, then at 
the lowest allowed antineutrino energies, only the for- 
ward positrons will survive the cut, biasing (cos 8) to 
be larger. (A similar effect occurs due to the energy- 
dependent efficiency for positron detection.) These and 
related effects increase the expected value of (cos 8) to 
0.16, in good agreement with observation f30| . This sug- 
gests that LSND is indeed observing D e + p — > e+ + n 
events, whatever the origin of the P e . 



E. Relation to neutron beta decay parameters 

It is worth noting that the parameter a(E„) is the cor- 
relation coefficient a between the positron and antineu- 
trino momenta in neutron beta decay pi] ]. The pa- 
rameter a is difficult to measure in neutron beta decay, 
since the antineutrino momentum can be reconstructed 
only by measuring the proton recoil momentum. The 
last measurement [[52|, more than 20 years ago, yielded 
a = —0.102 ± 0.005. It is tempting to speculate that it 
might be alternatively measured with D e + p — > e + + n 
via measurement of the positron energy and angle. 

It must be pointed out, however, that even at mod- 
est E v , weak magnetism is an important correction, and 
thus a measurement of a(E v ) would probe a combina- 
tion of /, g, and /2, thus providing instead a possible 
test of weak magnetism. To pursue this speculation fur- 
ther, consider an experiment in which (cos 8) could be 
measured at a fixed E v , and consider how well the value 
°f (/ + fs) could be measured. That is, an experimental 
test of whether the value of (/ + /2) extracted matches 



the value (4.706) predicted by the Conserved Vector Cur- 
rent hypothesis |3^]. The expectation for (cos 8) is ap- 
proximately given by Eq. (|23|). From statistics alone, 
5 ({cos 9)) ~ 1/V3N for N events, so that 



S(f + h) 



d(f + h) 
d (cos 8) 



8((cos8)) 



1 



E v /MyfN 



(27) 



We assume the standard V — A theory and neglect 
the (small) form factor variation. Some previous tests 
of weak magnetism were made by measuring extremely 
small distortions in the beta spectra of the A = 12 and 



A = 20 systems (see Refs. |34,|35| for a review). These 
experiments reached a precision on (/-I-/2) of about 10%. 

As noted above, for a Galactic supernova at 10 kpc, 
observed in SuperKamiokande, ~ 10 4 v e + p — > e + + n 
events are expected. The time distribution of events is of 
course irrelevant. The spectrum f(E v ) is also irrelevant 
insofar as it will be determined from the data, since the 
measured E e and cos6> e can be used to reconstruct E v 
for each event. At each value of E v , (cos 8) and hence 
(/ + fi) could be measured. While the typical energy 
expected is ~ 20 MeV, about 2/3 of the events are at 
higher energies. Thus one might plausibly expect that 
a test of weak magnetism at the ~ 20% level might be 
made (note also that the error varies linearly with the 
assumed supernova distance). 



IV. THE NEUTRON ANGULAR DISTRIBUTION 
AND APPLICATIONS 

It is also of interest to consider the angular distribu- 
tion of the neutrons from v e + p — > e + + n, since the 
neutrons are often detected as well. In fact, the obser- 
vation of the neutron capture in a delayed coincidence 
with the positron is the distinguishing signature of the 
antineutrino interaction with protons, allowing suppres- 
sion of backgrounds. 

There is an angular correlation between the antineu- 
trino direction and the initial direction of the neutron. 
Since in the laboratory system the proton is at rest, mo- 
mentum conservation requires that 



Pv = Pe 



Pn 



Also, 



be I < V(E» ~ A) 2 



< E„ 



(28) 



(29) 



so that the neutron must always be emitted in the for- 
ward hemisphere. In fact, the maximum angle (8 n )max 
between the antineutrino and initial neutron directions 
is achieved when the neutron and and positron momenta 
are perpendicular. Neglecting 0(1/ M), this is 



j2E v & - (A 2 



/ n J max 



E„ 



(30) 
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FIG. 4. The average neutron cosine and the cosine of the 
maximum neutron angle versus antineutrino energy. Note 
that cos(8 n ) max — » 1 both at threshold and as E v — > oo. 



At threshold, the neutron direction is purely forward, and 
at reactor energies, still largely so. In Fig. 4 we plot the 
quantity cos(9 n ) max as well as the average (cos(0 n )), both 
evaluated numerically, where the latter was weighted 
with the differential cross section. At 0(1/M) (see 
Eq. (Il3])), the neutron kinetic energy is 



T n = 



E V E, 



(0) 



(31) 



It is often possible to localize, at least crudely, the 
points where positron annihilated (essentially the point 
of creation for targets with an appreciable density) and 
where the neutron was captured. Even though the neu- 
tron is captured only after many elastic scatterings, its 
final position maintains some memory of its initial di- 
rection, as we now show. For a monoenergetic source of 
neutrons moving initially along the x-axis, the distribu- 
tion of final positions is Gaussian distributed, with equal 
widths a x , a y , and a z . The average final position (x) is 
displaced from the origin, and (y) = (z) = 0. 

We show the results of our Monte Carlo simulation 
in Fig. 5, implemented by following the principles given 
by Fermi |36j. We assumed a liquid of (CH2) n , with or 
without Gd doping (0.1% by mass), and a density of 0.80 
g/cm 3 . The results in Fig. 5 can easily be rescaled to an- 
other density by multiplying by (0.80/p). Neutrons are 
moderated by elastic scattering until they reach thermal 
energy. At thermal energy, elastic scattering changes the 
neutron direction, but, on average, not its energy. We 
also implement capture on protons and Gd. The calcu- 
lated capture times on undoped and doped scintillator 
are in good agreement with expectation. 

The most significant input is the fact that the neutron 
elastic scattering cross sections are almost independent 
of energy from about 10 keV down to about 0.1 eV |p7|. 
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FIG. 5. The shift (solid line) (x) and width (long-dashed 
line) a x — yl (x 2 ) — (x) 2 for monoenergetic neutrons (initial 
kinetic energy T n ) emitted from the origin, moving initially 
along the :r-axis. Note (y) — (z) = 0, and a x — a y = o z . We 
used a (CH2) n liquid of density 0.80 g/cm 3 , with or without 
0.1% Gd doping by mass. 



For the first ~ 5 scatterings, the neutron maintains some 
sense of its initial direction, with the angular distribu- 
tion at a given scattering depending only on the ratio of 
incoming and outgoing energies. This stage determines 
(x) ~ 1.7 cm, the average displacement of the final point 
from the starting point. All subsequent scatters only en- 
large the size of the neutron cloud. The width a x does de- 
pend on the the initial neutron energy because at higher 
energies more scatterings are required to moderate the 
neutron to thermal energy. 

Above 10 keV, the variation of the cross sections with 
energy should be taken into account; doing so would in- 
crease (x) and a x somewhat for T n ~ 100 keV and more 
substantially for higher energies. Below about 0.1 eV, 
there is a chemical binding effect depending on the mod- 
erating material that increases the cross section; taking 
that into account would make a x smaller. Finally, from 
Fig. 4, the struck neutron is not exactly forward, as as- 
sumed, but has cos 8 n ~ 0.9 for reactor energies. Taking 
that into account would reduce our shift to (x) ~ 1.5 cm, 
in agreement with the results of Ref . . 

In fact, in the Gosgen reactor antineutrino experi- 
ment [3S| l the neutron displacement was clearly observed, 
at ~ lOcr level p0| . This was possible because the detec- 
tor was composed of alternating walls of scintillator and 
3 He neutron detectors. For a given wall of scintillator 
in which the reaction occurred, and the positron was de- 
tected, more neutrons were observed in the 3 He slab away 
from the reactor than towards the reactor (in fact, the 
ratio was about 2:1). A similar effect was observed |4jJ| 
in the Bugey 3 experiment p2[, also using a segmented 
detector. 

The neutron-positron separation is also being used |3q] 
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by the Chooz experiment [J43J], which is based on an un- 
segmented detector. The neutron position is only de- 
tected with a precision of about 20 cm, but nevertheless a 
statistically significant displacement of positron and neu- 
tron detection positions along the antincutrino direction 
is seen. 

Given a reliable calculation of the neutron transport in 
the detector, and hence the expected neutron distribu- 
tions, this technique would allow a direct determination 
of the detector background from the measured asymme- 
try. Such an analysis is being pursued for the Palo Verde 
reactor experiment, and a forward-backward asymmetry 
between different cells is seen in the current data E3| . 



V. CONCLUSIONS 
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We have given an expression for the differential cross 
section in the positron cos 9, valid to order 1/M. Recoil 
and weak magnetism corrections have a large energy- 
dependent effect on the positron angular distribution, 
changing it from slightly backward at low energies to 
isotropic at about 15 MeV and slightly forward at higher 
energies. 

Our main result for the total cross section, valid to 
order 1/M, is obtained by integrating Eq. (|l4|). At low 
energies, this agrees with the well-known M — > oo result. 
At high energies, where the threshold can be neglected, 
this is in good agreement with Eq. (3.18) of Ref. B, which 
contains all orders in 1/M but assumes A — 0. Using 
our result, we have determined the largest A-dependent 
terms missing in the results of Ref. Q . The most accu- 
rate formula for the total cross section at all energies is 
obtained by using the result of Ref. || with our modifica- 
tions. The modified formula is essentially identical to our 
main result for E v < 30 MeV and is in good agreement 
with it (and the unmodified result of Ref. Q ) at higher- 
energies. 

The positron angular distribution is well-described by 
(cos#). Our main result, Eq. (p2]), valid to order 1/M, is 
an excellent approximation over the entire energy range 
considered (and in fact up to about E v ~ 150 MeV). 

A number of experimental applications are discussed 
in which the correct angular distribution is necessary for 
separating v e + p — > e + + n events from other reactions 
and from detector backgrounds. 

The neutron angular distribution is initially strongly 
forward. The random walk of the neutrons as they are 
moderated acts to erase this. However, the centroid 
of the final distribution of neutron capture positions is 
shifted in the direction of the initial motion. We discuss 
how this can be exploited experimentally. 



APPENDIX A: KINEMATIC RELATIONS 

In the center of momentum frame, the threshold is de- 
fined by the positron and neutron being produced at rest, 
so 



j^thr. 



(M n + m e ) 2 - Ml 
2(M„ + m e ) 



1.803 MeV. 



(Al) 



In the laboratory frame (where the proton is at rest), the 
threshold is 



E, 



thr. 



(M n +m e ) 2 -M* 

IMri 



1.806 MeV. 



(A2) 



Labeling the 4- momenta as v e (p u ) +p(p p ) —* e + (p e ) + 
n(p n ), we define the Mandclstam variables as 



s = (pv+Pp) 2 
t = (p u -p e ) 2 = M 

\2 7i,r2 i "1 



Ml 

f-2 



(A3) 



2M P E V 



M 2 p - 2M V {E V - E e ) (A4) 



(pu - Pn) = M; + mi 



2M p E e , 



(A5) 



evaluated in the laboratory frame, where we can also 
write t = q 2 = m\ — 2E v E e (l — v e cos8). 

The differential cross section in t can be written as 



da _ g| cos 2 C \M\ 2 

dt 7T 



Ml) 2 



(1 + A R ), 

V 1 inner / ' 



(A6) 



where \A4\ 2 is the amplitude squared (averaged over ini- 
tial spins, summed over final spins). This can be written 
as the differential cross section in the positron cos0 in 
the laboratory by using the Jacobian 



dt 



dcosO 



2E„pW 



M 



,(0) 



cos 8 



O 



1 

M 2 

(A7) 



The differential cross section in E e can be obtained by 
using dq 2 /dE e — 2M p . 
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